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Abstract. Many systems that can be described in terms of diffusion- limited 'chem- 
ical' reactions display non-equilibrium continuous transitions separating active from 
inactive, absorbing states, where stochastic fluctuations cease entirely. Their critical 
properties can be analyzed via a path-integral representation of the corresponding 
classical master equation, and the dynamical renormalization group. An overview 
over the ensuing universality classes in single-species processes is given, and gen- 
eralizations to reactions with multiple particle species are discussed as well. The 
generic case is represented by the processes A ^ A + A, and A — *• 0, which map onto 
Reggeon field theory with the critical exponents of directed percolation (DP). For 
branching and annihilating random walks (BARW) A — > (m + 1)A and A + A — > 0, 
the mean-field rate equation predicts an active state only. Yet BARW with odd m 
display a DP transition for d < 2. For even offspring number m, the particle number 
parity is conserved locally. Below d' c ~ 4/3, this leads to the emergence of an inac- 
tive phase that is characterized by the power laws of the pair annihilation process. 
The critical exponents at the transition are those of the 'parity-conserving' (PC) 
universality class. For local processes without memory, competing pair or triplet 
annihilation and fission reactions kA — > (k — I) A, kA — > (k + m)A with k = 2,3 
appear to yield the only other universality classes not described by mean-field the- 
ory. In these reactions, site occupation number restrictions play a crucial role. 

PACS: 64.60. Ak, 05.40.-a, 82.20.-w 

1 Introduction: Active to absorbing state transitions 

The characterization of non-equilibrium steady states constitutes one of the 
prevalent goals in present statistical mechanics. Unfortunately, away from 
thermal equilibrium one cannot in general derive even stationary macroscopic 
properties from an effective free energy function. One might hope, however, 
that such a classification in terms of symmetries and interactions becomes fea- 
sible near continuous phase transitions separating different non-equilibrium 
steady states: Drawing on the analogy with equilibrium critical points, one 
would expect certain features of non-equilibrium phase transitions to be uni- 
versal as well, i.e., independent of the detailed microscopic dynamical rules 
and even the initial conditions. The emerging power laws and scaling func- 
tions describing the long-wavelength, long-time limit should then hopefully 
be characterized by not too many distinct dynamic universality classes. Yet 
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studies of a variety of non-equilibrium processes have taught us that critical 
phenomena as well as generic scale invariance far from thermal equilibrium, 
where restrictive detailed-balance constraints do not apply, are often consid- 
erably richer than their equilibrium counterparts. Indeed, intuitions inferred 
from the latter have frequently turned out to be quite deceptive. 

But as in the investigation of equilibrium critical phenomena, field theory 
representations supplemented with the dynamical version of the renormal- 
ization group (RG) provide powerful methods to extract and systematically 
classify universal properties at continuous non-equilibrium phase transitions. 
Additional indispensable quantitative tools are of course Monte Carlo simu- 
lations, and other numerical approaches and exact solutions when available 
(the latter are usually restricted to specific one-dimensional systems). An- 
alytical and numerical methods usually supplement each other: Simulation 
results often call for deeper understanding of the underlying processes, but 
also rely on some theoretical background as a basis for data analysis. 

A prominent class of genuine non-equilibrium phase transitions separates 
'active' from 'inactive, absorbing' stationary states where, owing to the ab- 
sence of any agents, stochastic fluctuations cease entirely |1I2| . These occur 
in a variety of systems in nature; e.g., in chemical reactions which involve 
an inert state that does not release the reactants A anymore. Another 
example are models for stochastic population dynamics, combining, say, dif- 
fusive migration of a species A with asexual reproduction A — > 2A (with rate 
a), spontaneous death A — » (at rate /i), and lethal competition 2 A — ► A 
(with rate A). In the inactive state, where no population members A are left, 
clearly all processes terminate. Similar effective dynamics may be used to 
model certain non-equilibrium physical systems. Consider for instance the 
domain wall kinetics in Ising chains with competing Glauber and Kawasaki 
dynamics [3j. Here, spin flips TTII^TTTI an d TTIT ^ TTTT ma Y be viewed 
as domain wall (A) hopping and pair annihilation 2 A — > 0, whereas spin ex- 
change f T4I —> T4TI represents a branching process A — > 3A. Notice that the 
para- and ferromagnetic phases respectively map onto the active and inac- 
tive 'particle' states, the latter rendered absorbing if the spin flip rates are 
computed at zero temperature, thus forbidding any energy increase. 

The simplest mathematical description for such processes uses kinetic rate 
equations, which govern the time evolution of the mean 'particle' density n(t). 
For example, the above population model leads to Fisher's rate equation 



It yields both inactive and active phases: For a < fi we have n(t — > oo) — > 0, 
whereas for a > // the density eventually saturates at n s = (a — fi)/X. The 
explicit solution (with initial particle density no) 



shows that both stationary states are approached exponentially in time. The 
two phases are separated by a continuous non-equilibrium phase transition at 



n(t) = {a-n) n(t) - A n{tf . 



(1) 




(2) 
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a = n, where the temporal decay becomes algebraic, n(t) = no/(l+noXt) — ► 
l/(Ai) as t — > oo, independent of the initial density. But Eq. (JTJ represents 
a mean-field approximation, as we have in fact replaced the joint probability 
of finding two particles at the same position with the square of the mean 
density. As in equilibrium, however, critical fluctuations are expected to in- 
validate simple mean-field theory in sufficiently low dimensions d < d c , which 
defines the upper critical dimension. A more satisfactory treatment therefore 
necessitates a systematic incorporation of spatio-temporal fluctuations, in- 
cluding specifically the particle correlations as induced by the dynamics. 



2 From the master equation to stochastic field theory 

The renormalization group study of (near-) equilibrium dynamical critical phe- 
nomena relies on phenomenological Langevin equations for the order param- 
eter and 'slow' hydrodynamic variables associated with conservation laws 
0j . All other degrees of freedom are treated as Gaussian white noise, whose 
second moment is related to the relaxation coefficients through Einstein's 
fluctuation-dissipation theorem. As we shall see, however, such a description 
is not necessarily possible at all in reaction-diffusion systems. To the very least 
one would have to invoke fundamental conjectures on the noise correlators; 
but far from equilibrium these often crucially influence long-wavelength prop- 
erties. It is therefore desirable to construct a long-wavelength effective theory 
for stochastic processes directly from their microscopic definition, without re- 
course to any serious additional assumptions or approximations. 

Fortunately, there exists an established procedure to derive the Liouville 
time evolution operator for locally interacting particle systems immediately 
from the classical master equation, wherefrom a field theory representation 
is readily obtained [Sj. The key point is that all possible configurations can 
be labeled by specifying the occupation numbers n% of, say, the sites of a d- 
dimensional lattice. Let us for now assume that there are no site occupation 
restrictions, i.e., any m — 0,1,2,... is allowed (we shall address effects of 
particle exclusions in Sec.EJl- The master equation governs the time evolution 
of the configurational probability P({rii}; t). For example, the corresponding 
contribution from the binary coagulation process 2A — ► A at site i reads 

d t P{n l -t)\ x = \\ i {n l + l)n l P{n t + l-t)-n l {n l -l)P{n l ]t)^ . (3) 

This sole dependence on the integer variables {rii} calls for a representation 
in terms of bosonic ladder operators with the standard commutation relations 
[o», a-] = Sij, and the empty state |0) such that Oj|0) = 0. We next define the 

Fock states via = ni( a !) ni |0) (notice the different normalization from 

standard quantum mechanics), and thence construct the formal state vector 



!*(*)> = £ p ({"*>;*) IKI) • (4) 
{»«} 
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The linear time evolution imposed by the master equation can be cast into an 
'imaginary-time Schrodinger' equation d t \$(t)) = —H \$(t)) with a generally 
non-Hermitian, local 'Hamiltonian' H ({aj}, {at})- For instance, the on-site 
coagulation reaction is encoded in this formalism via = — A (1 — aj) a\a 2 . 

Our goal really is to evaluate time-dependent statistical averages for 
observables F, necessarily also just functions of the occupation numbers, 
(F(t)} = J2{m} F(i n i}) P{{ n i\it)- Straightforward algebra using the iden- 
tity [e a , a^] = e a shows that this average can be written as a 'matrix element' 

(F(t)) = (V\ F({ ai }) = (V\ F({ ai }) e~ m |#(0)) (5) 

with the state vector \<P(t)) and the projector state (V\ = (0| FJi e ai ; notice 
that (V\0) = 1. Probability conservation implies 1 = (V\ e~ Ht |0(O)), i.e., for 
infinitesimal times (7 7 |^(0)) = 1 and (V\H — 0, which is satisfied provided 
H({a\ — * l},{ai}) = 0. We remark that commuting the factors e ai to the 
right has the effect of shifting all a\ — > 1 + a\. One may then emply coherent 
states, as familiar from quantum many-particle theory [S], to represent the 
matrix element (0 as a functional integral with statistical weight exp(— S). 
Omitting contributions from the initial state, the action becomes 

S[{$i}, m] = [dt[j2 i>i dti>i + H({^}, {^})] . (6) 
•* i 

Taking the continuum limit finally yields the desired field theory that faith- 
fully encodes all stochastic fluctuations. 

3 Reggeon field theory and directed percolation (DP) 

Let us now return to our population dynamics example with random walkers 
A (with diffusion constant D in the continuum limit), subject to the Gribov 
reactions A ^ A + A and A — > 0. The corresponding field theory reads 

S= Jd d xalt ^{d t -D\7 2 )^ + a{l-^)Hj-fi(l-4>)4>~X{l-^tlj 2 .(7) 

The stationarity condition SS/Sip = is always solved by ip = 1; upon 
inserting this into SS/dtjj — 0, and identifying n(t) = (tp(x,t)), one recovers 
Fisher's mean-field rate equation (Q. Higher moments of the field ip, however, 
cannot be immediately connected with density correlation functions. In terms 
of an arbitrary momentum scale k, we record the naive scaling dimensions 
[<0] = k°, [ip] — n d , [a] — k 2 — [//], and [A] = n 2 ~ d . Hence the decay and 
branching rates constitute relevant operators in the RG sense, whereas the 
annihilation process becomes marginal at al — 2. Next we expand the action 
about the stationary solution ip = 1, i.e., introduce ip(x,t) = %j)(x,t) — 1, 
whereupon we arrive at 

S = I al d x dt \<<p(dt - L>V 2 )^ + (m - o-)4>ip -ai> 2 ip + X ip(l + ^)^ 2 j • (8) 
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Inspection of the one-loop fluctuation corrections shows that the effective 
coupling is in fact u = V o~X, with scaling dimension [u] = n 2 ~ d / 2 , whence 
d c = 4. At least in the vicinity of d c , X becomes irrelevant; certainly the ratio 
X/u scales to zero under subsequent RG transformations. Simple rescaling 
ip = 4>yJ X/a, tp — (f>yj a jX then leaves us with the effective field theory 



S cff [<M = d d xdt 



0[a t + ^(r-V 2 )]0 + U (^ 2 ~0 2 ^)] , (9) 

where r = (/i — <j)/D. This should capture the critical behavior for the non- 
equilibrium phase transition at r = in our population dynamics model. 

The action @ is known in particle physics as ' Reg g eon' field theory [7]. It is 
invariant with respect to 'rapidity inversion' <p(x,t) — -» —<p(x,—t), 4>{x,t) — > 
—<f>(x, —t). Quite remarkably, the very same action is obtained for the thresh- 
old pair correlation function [H] in the geometric problem of directed perco- 
lation (DP) U3- Fig.^b) depicts a critical directed percolation cluster, con- 
trasted with the structure emerging at the threshold of ordinary isotropic 
percolation. At d c = 4, the critical exponents as predicted by mean-field the- 
ory acquire logarithmic corrections, and are shifted to different values by the 
infrared-singular fluctuations in d < 4 dimensions. By means of the standard 
perturbational loop expansion in terms of the diffusion propagator and the 
vertices oc u, and the application of the RG, the critical exponents can be 
computed systematically and in a controlled manner in a dimensional expan- 
sion with respect to e — 4 — d. The one-loop results, to first order in e, as well 
as reliable values from Monte Carlo simulations in one and two dimensions 
are listed in Table ^ Moreover, as a consequence of rapidity invariance 
there are only three independent scaling exponents, namely the anomalous 
field dimension n, the correlation length exponent v, and the dynamic critical 
exponent z. All other exponents are then fixed by scaling relations, such as 

f3 = ^(z + d-2 + r])=zva (10) 

for the order parameter and critical density decay exponents. 

It is worthwhile noting that Reggeon field theory can be viewed as a 
dynamic response functional and therefore is equivalent to an effective 
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Table 1. Critical exponents for the saturation density n s , correlation length £, 
time scale t c , and critical density decay n c (t) for the DP universality class 

DP exponents d = 1 d = 2 d = 4 — e 

n a ~ |r|" /?« 0.2765 /3 ~ 0.584 /3 = 1 - e/6 + 0(e 2 ) 

i/» 1.100 i/« 0.735 z/ = 1/2 + e/16 + 0(e 2 ) 

t c ~ £ z ~ |r|- z " z « 1.576 2 » 1.73 2 = 2 - e/12 + 0(e 2 ) 

n c (t)~t~ a a « 0.160 a « 0.46 a = 1 - e/4 + 0(e 2 ) 



Langevin equation. To this end, we integrate out the field 0, which yields the 
statistical weight exp(— G) with 

G[<f>] = Jd d xdt d t (t> + D(r~V 2 )(j) + u<j) 2 ^Au(p. (11) 

After rescaling, we may interpret (|ll(l as the Onsager-Machlup functional 
associated with the Gaussian noise distribution for the stochastic process 

d t n(x, t) = D(V 2 - r) n(x, t) - A n(x, tf + £(x, t) , (12a) 
(C(x,t))=0, (((x,t)((x',t')) = 2an(x 1 t)6(x-x')5(t-t') . (12b) 

Here, the absorbing nature of the inactive state is reflected in the fact that 
the noise correlator is proportional to n(x,t). Of course, Eq. (|12b(l really 
means that the local density is to be factored in when the noise average is 
taken, c.f. Eq. (| 1 1 |t . Alternatively, One may define ^(a; , t) = y/n(x, t) rj(x, t), 
whereupon (rj(x,t)r)(x',t')) = 2aS(x— x') S(t— t'), at the cost of introducing 
'square-root' multiplicative noise into the Langevin equation l|12afl . Within 
the Langevin framework, we can readily generalize to arbitrary reaction and 
noise functionals r[n] and c[n]: 

d t n(x, t) = DV 2 n(x, t) - r[n(x, t)] + £(x, t) , (13a) 
(C(x, t))=0, (C(x, t) C(x', t')) = c[n{x, t)] S(x - x 1 ) S(t - t') . (13b) 

In the spirit of Landau theory, we may then expand the functionals r[n] and 
c[n] near the inactive phase (n -C 1). In the absence of spontaneous parti- 
cle production, both must vanish at n = 0, which is the condition for an 
absorbing state. Keeping only the lowest-order, relevant terms in the expan- 
sions with respect to n, we thereby infer that any active to absorbing state 
phase transition in a single-species system should genetically be described 
by Eqs. (|12afl and (|12b|l . i.e., Reggeon field theory ©. Consequently, in the 
absence of any special symmetries, memory effects, and quenched disorder, 
we expect to find the critical exponents of the DP universality class ^21 ■ 
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Fig. 2. Mean-field phase diagram for the two-stage unidirectionally coupled DP 
process. The arrows indicate active to absorbing transitions for the A and B species. 
The dotted parabola marks the boundary of the multi-critical regime |14| 

4 Variants of directed percolation processes 

In fact, DP-type processes with even arbitrarily many particle species have 
been fully classified. Consider the reactions A ^ A + A, A ^ 0, B x± 
B + B, B — ► 0, etc., supplemented with bilinear couplings A — > B + B, 
A + A — > B, . . . Higher-order reactions then turn out to be irrelevant in 
the RG sense, and remarkably the critical behavior of the multi-species DP 
system is once again governed by the DP fixed point ^21- However, if we allow 
for unidirectional linear couplings through particle transmutations A —> B, 
B — ► C, . . ., with rates fi AB , P^BCi ■ ■ ■> multi- critical behavior may ensue |14|. 
This becomes already manifest on the level of the coupled rate equations 

d t n A = D (V 2 - r A ) n A - \ A n\ , (14a) 
d t n B = D (V 2 - r B ) n B — X B n 2 B + [i AB n A , (14b) 

etc. For as long as the A species is in the active phase (r A < 0), the B 
particle density will be non-zero as well. As depicted in Fig. El this effectively 
'folds' half of the decoupled B transition line (r A < 0, r B = 0) over onto 
r A = 0, r B > 0. Along this half-line, the B particles become 'enslaved' by 
the A species; and so forth further down the hierarchy of particle species. As 
r A — > and r B — * simultaneously, one encounters a non-equilibrium multi- 
critical point. While the DP exponents v and z governing the correlation 
length and critical slowing down remain unchanged, one finds successively 
reduced values for the order parameter exponents on the jth hierarchy 
level 1 as listed in Table |3 The crossover exponent associated with the 
multi-critical point is <p = 1 to all orders in the e expansion 13 . 

Another mechanism to induce a different universality class in a two-species 
system is to link diffusing agents A with a passive, spatially fixed, and initially 
homogeneously distributed species X through the reactions X + A — > A + A 
and A — ► 0. One may then integrate out the X fluctuations; upon expanding 
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Table 2. Simulation and one-loop RG results for the saturation density critical ex- 
ponents on the first three hierarchy levels in unidirectionally coupled DP processes 





d = 1 


d= 2 


d = 3 


d = 


4 - e 




0.280(5) 


0.57(2) 


0.80(4) 


1 - 


e/6 + 0(e 2 ) 




0.132(15) 


0.32(3) 


0.40(3) 


1/2 


-e/8 + 0(e 2 ) 


0(3) 


0.045(10) 


0.15(3) 


0.17(2) 


1/4 


-O(e) 



about the mean-field solution, the resulting effective action becomes [To\ 



SefffcM = hrxdt 



(x,t) [d t +D(r-V 2 )] </>(x,t) 

+2Du4>(x,t)(/>(x,t) J 4>(x,t')dt' -ucf>{x,t) 2 <p(x,t) . (15) 

The elimination of the passive particles X thus induces memory of all preced- 
ing times in the particle annihilation vertex. The DP rapidity inversion in- 
variance is replaced by its non-local counterpart 4>(x, t) — > —D^ 1 dt<j>(x, —t), 
4>(x,t) -> —D J (f>(x,t')dt'. The upper critical dimension for this dynamic 
percolation universality class is shifted to d c — 6. Its static critical exponents 
are precisely the ones that characterize a critical isotropic percolation cluster 
|15| . compare Fig. Qfa). To first order in e = 6 — d, one finds rj = — e/21, 
v = 1/2 + 5e/84, z = 2 - e/6, and /3 = i/(d - 2 + r?)/2 = 1 - e/7. Multi- 
species generalizations proceed in the same way as for DP, with similar re- 
sults: Whereas non-linear couplings to other particle species preserve the 
dynamic universality class, a multi-critical point emerges for unidirectional 
particle transmutations, with crossover exponent (p — 1 13 . 



5 Diffusion-limited annihilation processes 



As a preparation for the following Scc.[fjl let us now investigate the fcth order 
annihilation reaction kA — * 0. The associated rate equation reads dt n(t) = 
—\n(t) k . For radioactive decay (k = 1), it is naturally solved by the familiar 
exponential n(t) — uq e _At , whereas one obtains power laws for k > 2, namely 



n(t) = [nl~ k + (k-l)Xt] 



-i/O-i) 



(16) 



In order to consistently include fluctuations in the latter case, we again start 
out from the master equation, wherefrom we derive the action |16j 



S[i>, rl>]= d d x dt U (d t - DV 2 ) — X(l — *jj k ) ^ 



(17) 



After performing the shift i/j(x, t) = 1 + ip(x, t), it becomes evident that this 
field theory does not have a simple Langevin representation. For in order 
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to interpret ip as the corresponding noise auxiliary field, it should appear 
quadratically in the action only, and with negative prefactor. Thus even for 
the pair annihilation process, the Langevin equation derived from the action 
fTTjl would entail unphysical 'imaginary' noise with c[n] = —2 An 2 . 

The scaling dimension of the annihilation vertex is [A&] = k 2- ^ -1 ^, 
whence we infer the upper critical dimension d c {k) = 2/(k— 1). This leaves the 
possibility of non-trivial scaling behavior in low physical dimensions only for 
the pair (k = 2) and triplet (k = 3) processes. Analyzing the field theory (| 1 T|l 
further, we see that the diffusion propagator does not become renormalized 
at all. Consequently, r\ = and z = 2 to all orders in the perturbation expan- 
sion. The simple structure of the action permits summing the entire vertex 
renormalization perturbation series by means of a Bethe-Salpeter equation; 
in Fourier space it reduces to a geometric series of the one- loop diagram [TB). 

For pair annihilation, this yields the following asymptotic behavior for the 
particle density: n(t) ~ (At) , i.e., the reaction- limited power law of the rate 
equation for d > 2; but diffusion-limited decay n(t) ~ (Dt) -d / 2 for d < 2. 
At d c = 2, one finds the logarithmic correction n(t) ~ (Z?t) -1 lnDt. The 
slower decay for d < 2 originates in the fast mutual annihilation of any close- 
by reactants; after some time has elapsed, only well-separated particles are 
left. The annihilation dynamics thus produces cmti-correlations, mimicking 
an effective repulsive interaction (which actually provides the interpretation 
for the negative correlator c[n]). In the ensuing diffusion-limited regime, the 
typical particle separation scales as £(t) ~ (Z?t) -1 / 2 , whereupon indeed n(t) oc 
£(t)~ d ~ (Dt)~ d / 2 . The same power laws hold for the pair coagulation process 
2A — > A, albeit with different amplitudes. Replacing ordinary diffusion with 
long-range Levy flights with probability oc r -d-p of hopping a distance r 
(p < 2) results in n(t) ~ (Dt)~ d / p for d < d c = p ^7). For triplet annihilation, 
one can similarly show that the density decays as n(t) ~ (At) -1 / 2 for d > 1, 

1 /2 

with mere logarithmic corrections n(t) ~ [(Dt)^ 1 InDt] at d c = 1 |lf>| . 

Generalizations of the pair annihilation reaction to multiple particle types 
introduce interesting new physics. For the two-species case A + B — * 
(with no concurrent reactions of identical particles), the rate equations read 
dtnA/B = —Xuatib- With equal initial densities hao = n-so they are again 
solved by fT-A/s(t) ~ (At) -1 ; however, with hao > n,Bo, say, one obtains 
ns(t) ~ exp[— (hao — nso) At] for the minority species, while the majority 
density saturates at jia s > 0. In order to establish the effects of spatial fluc- 
tuations, it is crucial to notice that the density difference ua — tib remains 
strictly conserved under the reactions; for Da = Db it simply obeys the 
diffusion equation [IS) . Consequently, regions with A or B particle excess 
become amplified in time. As a result, when hao = ubo, one finds that for 
dimensions d < 4 species segregation into A/B rich domains occurs The 
annihilation processes are then confined to sharp reaction fronts, leading to 
a decelerated density decay n^/s(t) ~ (Dt) -d / 4 . For unbalanced initial con- 
ditions, stretched exponential relaxation ensues for d < 2: lnns(t) ~ — t d / 2 , 
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whereas lnns(t) ~ — tj lni at d c = 2 20 . In one dimension, special ini- 
tial configurations may change this picture: Consider, e.g., the alternating 
arrangement . . . ABABAB ... of particles that upon encounter react with 
probability one. Now there is no reason anymore to distinguish between A 
and B, and the system is in the 2 A — » universality class. 

An obvious question is then what happens for diffusion-limited pair an- 
nihilation of q > 2 particle species, with equal initial densities as well as 
reaction and diffusion rates In contrast with the two-species case, there 
exists no local conservation law. Furthermore the renormalization of the re- 
action vertex proceeds exactly as for 2 A — > 0. Consequently, at least for 
d > 2, where the initial state is not that crucial, the long-time limit should 
in fact generically be governed by the single-species pair annihilation univer- 
sality class |22] • This is obvious for q — > oo: In this limit, the probability of 
like particles ever meeting vanishes, which renders the distinction of differ- 
ent species meaningless. However, in one dimension, at least in the limit of 
large reaction rates (which should describe the asymptotic regime) , particles 
of different types cannot pass each other. This topological constraint allows 
for species segregation to occur. Indeed, a simplified deterministic version of 
the g-species pair annihilation process yields [22j 

n(t) ~ , with a(q) = (q - l)/2q , (18) 

which correctly reproduces a(2) = 1/4 and a(oo) = 1/2 in d = 1. The asymp- 
totic decay (|18|l along with the subleading correction ~ t -1 / 2 of the pair anni- 
hilation process without segregation were recently confirmed in extensive sim- 
ulations |2H| . Yet again, special initial conditions such as . . . ABC D ABC D . . . 
may prevent segregation and instead lead to the 2A — > decay law. 



6 Branching and annihilating random walks (BARW) 

In order to allow again for a genuine phase transition, we combine the an- 
nihilation kA — > (k > 2) with branching processes A — > (m + 1) A, The 
associated rate equation for these branching and annihilating random walks 
(BARW) reads d t n(t) = —\n(t) k + an(t), with the solution 

n(t) = n s /(l + [(n s /n ) k - 1 - l] e -(*-i)"*) 1/( * _1) . (19) 

Mean-field theory thus predicts the density to approach the saturation value 
n s = (er/A) 1 /^ -1 ) as t — > oo for any positive branching rate a. Above the 
critical dimension d c {k) — 2j(k— 1) therefore, the system only has an active 
phase; a c = represents a degenerate 'critical' point, with scaling exponents 
essentially determined by the pure annihilation model: a = l/(k — 1) = j3, 
v = 1/2, and z = 2. However, Monte Carlo simulations revealed a much 
richer picture, in low dimensions clearly distinguishing between the cases of 
odd and even number of offspring m |3I24| : For k = 2, d < 2, and m odd, 
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a transition to an inactive, absorbing phase is found, characterized by the 
DP critical exponents. On the other hand, for even offspring number there 
emerges a phase transition in one dimension, described by a novel universality 
class with a » 0.27, w 0.92, v ps 1.6, and z ps 1.75. 

The above mapping to a stochastic field theory, combined with RG meth- 
ods, elucidates the physics behind those remarkable findings j^S] • The action 
for the most interesting pair annihilation case becomes 



drxdt 



ip(dt—D V ) -0 - A (1 - ip ) ip +a(l — ip m ) ip ip , (20) 



which in general allows no direct Langevin representation. Upon combining 
the reactions A — > (m + 1)A and 2A — > 0, one notices immediately that the 
loop diagrams generate the lower-order branching processes A — > (m — 1)A, 
A — > (to — 3) A . . . Moreover, the one-loop RG eigenvalue = 2 — to(to + 1)/2 
(computed at the annihilation fixed point) shows that the reactions with 
smallest to are the most relevant. For odd to, we see that the generic situation 
is represented by to = 1, i.e., A — > 2 A, supplemented with the spontaneous 
decay A — ► 0. After a first coarse-graining step, this latter process (with rate 
n) must be included in the effective model, which hence becomes identical 
with the action Q . Thence we are led to Reggeon field theory © describing 
the DP universality class, provided the induced decay processes are sufficiently 
strong to render a c > 0. Yet for d > 2 the renormalized mass term ctr — fin 
remains positive, which leaves us with merely the active phase. For d < 2, 
however, the involved fluctuation integrals are infrared-divergent, thus indeed 
allowing the induced decays to overcome the branching processes to produce a 
non-trivial phase transition. As function of dimension, the critical exponents 
display an unusual discontinuity at d c = 2, as they jump from their DP to 
the mean-field values as a result of the vanishing critical branching rate |25| . 

ft is now obvious why the case of even offspring number to is funda- 
mentally different: Here, the most relevant branching process is A — ► 3A, 
and spontaneous particle death with associated exponential decay is not gen- 
erated, which in turn precludes the previous mechanism for producing an 
inactive phase with exponential decay. This important distinction from the 
odd-m case can be traced to a microscopic local conservation law, for the re- 
actions 2A — > and A — > 3A, A — > 5^4 . . . always destroy or produce an even 
number of reactants, preserving the particle number parity. Formally, this is 
reflected in the invariance of the action (|20|l under the combined inversions 
ip — > —ip, ip — > —ip. As we saw earlier, the branching rate a certainly consti- 
tutes a relevant variable near d c = 2. Therefore the phase transition can only 
occur at a c — 0, and for any a > there exists only an active phase, described 
by mean-field theory. In two dimensions one readily computes the following 
logarithmic corrections: £(er) ~ c -1 / 2 ln(l/cr), and n(a) ~ a [ln(l/cr)]~ 2 |23] , 

However, setting to — 2 in the one-loop value for the RG eigenvalue 
y a , we notice that fluctuations drive the branching vertex irrelevant in low 
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INACTIVE 
PHASE 



MEAN 
FIELD 



Fig. 3. Stationary states and unstable RG fixed point 1/g* for BARW with even 
offspring number (PC universality class) as function of dimension d |25| 



dimensions d < d' c « 4/3. More information can be gained through a one- 
loop analysis at fixed dimension, albeit uncontrolled (23 ■ The ensuing RG flow 
equations for the renormalized branching rate <jr = a/Dn 2 , and annihilation 
rate X R = C d \/Dn 2 - d , with C d = r{2 - d/2)/2 d ~ 1 ir d / 2 read (for m = 2): 



d(J R 

~dl= aR 



(1 + <J R ) 2 - d / 2 



-df = XR 



2-d- 



A, 



(1 + a R ) 2 - d l 2 



(21) 



The effective coupling is then identified as g — Xr/(1 + OR) 2 ~ d l 2 1 which 
approaches the annihilation fixed point g* = 2 — d as a R — > 0, whereas 
for or — > oo the flow tends towards the active state Gaussian fixed point 
<?0 = 0. The separatrix between the two phases is given by the unstable 
RG fixed point g* = 4/(10 — 3d), which enters the physical regime below the 
borderline dimension d' c ss 4/3, as shown in Fig. [31 For d < d' c , this describes a 
dynamic phase transition with a c > 0. The aforementioned fixed-dimension 
RG analysis yields the rather crude values v 3/(10 — 3d), z f» 2, and 
/3 w 4/ (10— 3d) for this parity-conserving (PC) universality class. The absence 
of any mean-field counterpart for this transition precludes a direct derivation 
of the 'hyperscaling' relations HlOjl. Amazingly, in this non-equilibrium system 
fluctuations generate rather than destroy an ordered phase (translating back 
from the domain wall to the spin picture) in low dimensions. The inactive 
state is characterized by a vanishing branching rate, and consequently by 
the algebraic pair annihilation density decay. For particles undergoing Levy 
flights, the existence of the power-law inactive phase is controlled by the 
anomalous diffusion exponent p, emerging for p > p c sa 3/2 in d = 1 26 . 

Invoking similar arguments for the case of triplet annihilation 3A — > 
combined with branching A — > (m + 1)A one would expect DP behavior 
with <j c > for m mod 3 = 1,2, as then the processes A — > 0, A — * 2 A, 
and 2A — * A are dynamically generated. For m = 3,6,..., on the other 
hand, there can be different, novel scaling behavior, but because of d c = 1 
it will be limited to merely logarithmic corrections in one dimension |25j . It 
is also interesting to generalize the even-offspring BARW to q species, such 
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that only equal particles can annihilate, Ai + Ai — » 0, but both reactions 
Ai — ► 3Ai (with rate a) and Ai — ► A4 + 2Aj (for j'^i, and with rate a') are 
possible. It turns out that the latter process always dominates, and in fact 
the ratio a /a' — ► under renormalization. Thus asymptotically one reaches 
the exactly analyzable q — > 00 limit, with a mere degenerate phase transition 
at branching rate a' c = 0. Below d c = 2, one finds the critical exponents 
a = d/2, = 1, v = 1/d, and z = 2 The situation for q = 1 is thus 
qualitatively different from all multi-component cases. 



7 Annihilation— fission reactions 

For single-species reactions without memory and disorder, the only remain- 
ing processes with potentially non-mean-field scaling behavior appear to be 
the combination of purely binary annihilation and fission reactions 2A — » A 
(with rate A) and 2A — > (m + 2) A (rate a m ) with d c = 2, and its triplet 
counterpart (d c = 1). The former reactions subsequently generate 2 A — > 
(m + 1)A, mA, . . ., 2(m + . . ., thus producing infinitely many couplings 
with identical scaling dimensions [22] ■ Upon including all these binary par- 
ticle production reactions, the phase transition is readily seen to occur at 
A c = J2m ma m- The inactive, absorbing phase (A > A c ) is obviously charac- 
terized by the power laws of the pure coagulation model. Yet for A < A c , the 
particle density diverges after a finite time, when no constraints on the site oc- 
cupation numbers n, are imposed. Thus the asymptotic density is finite only 
at the phase transition itself. These singular features of the 'bosonic' model 
with its highly discontinuous phase transition are overcome by restricting 
the site occupation numbers to m = 0, 1. Extensive Monte Carlo simulations 
have in fact revealed that this leads to a continuous transition, with criti- 
cal exponents that seem to belong to a novel universality class (pair contact 
process with diffusion, PCPD) with critical dimension d c = 2 28 . However, 
owing to the difficulty of obtaining truly asymptotic properties in this system, 
where reactions become extremely rare at low densities, the precise nature of 
this critical point in purely binary reactions has remained elusive and rather 
controversial. This applies even to density matrix RG studies [29] . 

It is therefore fortunate that recent work has demonstrated how to con- 
sistently implement site occupation restrictions into the bosonic field theory 
|.30| . For the above binary processes, the reaction part of the action becomes 

S = { d d xdt\a m (l-ip m ) V> 2 *P 2 e -(™+ 2 )^-A (!-</;) tp V 2 e- 2v ^} . (22) 



Here the exponential terms capture the occupation number limitations, with 
[v] = n~ d , which suggests that v represents a dangerously irrelevant coupling. 
Indeed, consider more generally the coupled reactions kA — > (k — I) A with 
< / < k and nA — > (n + m)A with n, m > 0, which display a continuous 
transition for k < n. The mean-field equations obtained from the associated 
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actions show that site occupation restrictions can be neglected at low densi- 
ties, yielding the critical exponents (3 = l/(fc — n), v = (k — l)/2(k— n), z = 2, 
and a = l/(fc — 1), except for the degenerate case k = n, where one finds 
mn s — ln(ma m /lX), whence j3 = 1, v = fc/2, z = 2, and a = 1/k [M132\ . For 
k = 1, expanding the exponentials leads to the action (JHJ, which establishes 
that the competing processes A — > 0, A — » 2/1, . . . with site exclusion yield a 
DP phase transition. The above field theory should also permit a systematic 
analysis of the fluctuation corrections for the purely binary and triplet reac- 
tions. For the latter, one expects mere logarithmic corrections at d c = 1 to 
the mean- field scaling laws; yet current simulations are inconclusive [32] . 

8 Concluding remarks 

In this overview, I have outlined how non-linear stochastic processes via their 
defining master equation can be represented by field theory actions, allowing 
for a thorough analysis and classification by means of the renormalization 
group. Systems with a single 'particle' species that display a non-equilibrium 
phase transition from an active to an inactive, absorbing state are generically 
captured by the directed percolation universality class. The second promi- 
nent example, sometimes applicable when additional symmetries (degenerate 
absorbing states) are present, is the parity-conserving universality class of 
even-offspring branching and annihilating random walks. The only other sce- 
narios for non-trivial critical scaling behavior appear to be provided by the 
solely pair or triplet annihilation-fission reactions, where site occupation re- 
strictions become relevant. The full classification of reaction-diffusion models 
with multiple particle species remains a formidable task. Even the difference 
in diffusivities may become a relevant control parameter . Specifically in 
one dimension, exclusion constraints can play a crucial role |34) . 

Another obviously important open problem concerns the influence of 
quenched disorder in the reaction rates. For example, a field theory RG inves- 
tigation for DP with random percolation threshold yields run-away flows |85| . 
reflected in intriguing simulation results with not entirely clear interpretation 
A very recent strong disorder RG study, supplemented with numerical 
density matrix RG calculations, has revealed a novel disorder fixed point |37|. 
A better understanding of spatially varying reaction rates might also explain 
the conspicuous rarity of clear-cut experimental realizations even for the sup- 
posedly ubiquitous DP universality class [SHI- In fact the single verification 
of DP scaling behavior appears to be its observation in spatio-temporal in- 
termittency in ferrofluidic spikes |39| . Thus many intriguing issues are still 
open; but I expect that in addition to increasingly more extensive Monte 
Carlo simulations and sophisticated numerical techniques, field theory rep- 
resentations and subsequent analysis by means of the renormalization group 
will remain an invaluable tool for the further understanding of cooperative 
phenomena and scale-invariance in interacting non-equilibrium systems. 
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